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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE No. 1712 


ANALYTICAL METHOD OP DE TE RMINING TRANSMISSION CP PARTICLES 
AND RADIATION THROUGH GREAT THICKNESSES CEF MATTER 
By A. R. Bobrovsky 


SUMMARY 

A matrix method of determining transmission, scattering, «nfl 
absorption of radiation, such as gamma rays ana neutrons, in large 
thicknesses of matter is presented. Experimental data as yet 
unavailable are required for complete utilization of the method. 
Material radioactivity whose intensity may be considered to be 
independent of time is also treated. 


INTRODUCTION 

An analytical method, developed at the NACA Cleveland labora- 
tory, is presented that can be used to determine the relative effi- 
ciency of a number of composite shields intended for use with one 
or more types of radiation. The procedure described is also useful 
in designing such shields. Application of the method is limited by 
the unavailability of data on the scattering and absorption of 
radiation by particular chemical elements in certain ranges of the 
energy spectrum. 

The matrix methods employed in handling a number of simul- 
taneous equations were used (reference 1) to solve a thermal- 
shielding problem analogous to the first case of this analysis. In 
reference 1, the efficiency of a number of thin shields, used with 
a source of thermal radiation, such that all radiation is either 
reflected at the surfaces of the shields or transmitted without 
absorption is determined; monochromatic radiation only is considered, 
and this radiation is assumed always normal to the shield. 

Shielding involving interactions between radiation and matter 
within the body of the shield, such as occurs with gamma rays and 
neutrons, is calculated herein. The following conditions, to which 
no experimental exceptions have been found for these types of 
radiation, were placed on the analysis: 
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(1) The amounts of -radiation transmitted, absorbed, and scattered 
in a given' shield are directly proportional to. the amount of radia- 
tion entering the shield for a given type of energy spectrum. Two 
energy spectrums are defined to he of the same type when the ratios 

of corresponding intensities at each energy are equal. This condi- 
tion may he alternatively stated as postulating zero saturation 
effect over the complete range of intensities and energies con- 
sidered. 

(2) The fractions of radiation of any energy that are transmitted, 
absorbed, and scattered in a given shield are independent of the 
presence of radiation of another energy. This condition enables the 
analysis to be made with quantities that obey laws of linear com- 
bination. 

(3) An elementary thickness of material cannot offer infinite 
resistance to the transmission of incident radiation in its original 
direction at its original energy. 

When these conditions were established, analysis proved feasible 
for cases in which polychromatic radiation was partly absorbed, 
partly degraded in energy, and partly transmitted unchanged in 
energy. Extensions of the concepts employed permitted consideration 
also of scattering at an angle to the direction of the incident beam, 
upgrading of energy as by nuclear reaction, and excitation of neutron 
emission by gamma-ray bombardment and of gamma rays by neutron 
bombardment. The scattered radiation from great thicknesses of 
material must be particularly considered not only because the mate- 
rial gives rise to much scattered radiation and secondary radiation, 
but because the shielding is required to cover a source that radiates 
spherically, so that the power of radiation at any point outside of 
the shielding is the integrated power from all points on the surface 
of the shield visible at the point of measurement. 

The method of calculation presented requires complete informa- 
tion on scattering and absorption in order to yield a true answer. 

If information on only the primary beam is available, the shielding 
calculation based on it must be in error by the secondary and 
scattered radiations that are being neglected. If the energy spec- 
trums of scattered radiation are unknown, the effect of energy on 
absorption must be neglected. Information is therefore desirable 
on the absorption and energy- Intensity distribution for radiation 
of all wavelengths to be encountered for all materials that will be 
handled in order to obtain the truest possible shielding analysis. 

If appreciable quantities of gamma radiation are converted to other 
radiation such as neutrons, these quantities must also be known. 
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A simplification in the data requirements is that data need, he 
available for only one elementary thickness of absorber and not a 
number of thicknesses. An elementary thickness is any given thick- 
en ness for •which absorption and scattering data are known or can 

ui conveniently be obtained, as, for example, 1 centimeter. Calcula- 

tions can be made by the method presented using currently available 
data, but the results obtained from the shielding analysis will be 
somewhat inaccurate. 

The data required for neutron absorption are those concerned, 
with inelastic scattering, elastic scattering, and the products of 
neutron-induced reactions. The data required here are similar to 
those, required for gamma rays eznept that they include not only the 
intensities of neutrons of different energise and dir ections but 
also the intensity and the direction of the energy spectrum of gamma 
radiation produced if appreciable. A number of parameters, such as 
grain size of metals, that have no effect on absorption of gamma 
radiation appear to influence neutron absorption (reference 2) and 
must be considered. 

Different materials are best for absorption at different 
incident spectrums of energy. Inasmuch as the spectrum varies 
through the shield, the shield of optimum absorption character- 
istics would contain different materials through its thickness. 

A suitable method of analysis consequently must be capable of 
determining the spectrum at a station within the shield, selecting 
an optimum absorber, and ascertaining the new emergent spectrum 
for the next station. 


S2MB0LS 

The following symbols are used in this method: 

Si i’kk element of shield 

. Ii power of radiation from (i-l) tlx element 

incident on i th element 

Ii + i power of radiation from i^ element 

incident on (i+l)’*' 11 element 

Ri power of radiation from i^ element 

incident on (i-l)^ element 


I 
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power of radiation from (i+l)^* 1 element 
incident on i^ element 

one-half of power of radiation from radio- 
activity in shield 

power of radiation from source of radia- 
tion to shield 

power- transmission' matrix for complete 
shield 

_ / 

power- transmission matrix for one element 
of shield 

matrix element 

power- transmission coefficient, or power- 
scattering coefficient 

power hack- scattering (reflection) 
coefficient 

power- conversion coefficient, radiation 
to mass and heat 

power-conversion coefficient, radiation 
to heat 

power-conversion coefficient, radiation 
to mass 

power of- radiation converted to mass and 
heat 

power of radiation converted to heat 

power of radiation converted to mass 

power-transmission, hack-scattering 
•(reflection), heat- conversion, mass- 
conversion, and total-conversion operators, 
respectively 

\ 

arbitrary single-valued function of inde- 
pendent variable 

arbitrary numerical coefficient 
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The following subscripts are used to designate positions anfl locations: 

a > &} • • •> v zones for radiation scattered at angles to 

primary direction 

1, 2, • • ., i, • • ., n successive elements of set (shielding or 

energies) 

The following superscripts are used to designate type of radiation: 
y gamma radiation 

n neutron radiation 

Y n conversion from gamma radiation to neutron 

radiation 

nY conversion from neutron radiation to 

gamma radiation 


"ANALYSIS 

The analysis of transmission of particles and radiation through 
matter may he conveniently "broken down into a number of separate 
cases. Each succeeding case contains an additional consideration 
of some phenomenon encountered in the transmission problem. 


Case I - Radioactive Shield 

The following hypothesis is made: 

Monochromatic radiation 
Radiation normal to shield 

Source of radiation both outside shield and from long- 
lived radioactivity of shield material 

Method. - The shield may be considered as consisting of elements 
orthogonal to the radiation. The i th element in a shield is Sj. 
Radiations from and S^ + p, nam ely and R^ + p, respectively. 



are Incident on S^. From emerge radiations I i+ 2 and 
to S i+ ^ and respectively. Badiation is assumed to origi- 

nate from sources outside the shield as veil as from an Independent 
source of radiation 2P^ within the shield. In practice, such an 

independent source might exist in the form of long-lived radio- 
isotopes in the shielding of a cyclotron that had operated over a 
period of time. 



Let t^, the power-transmission coefficient, he a number that 
determines the amount of power transmitted through the element; for 
a given thickness, it is a function of the wavelength of the radiation 
and of the material constituting the element. 



S imi la r ly the hack-scattering, or reflection, coefficient Is 



If the shield power 2P^ is distributed equally between the two 
faces of the element, then 

^d+l “ t i I i + r i R i+l + P I (3) 

R i * t i B i+l + r i x L + P I 


( 4 ) 
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If equation (3) is rewritten, 

h. - ^ *1+1 - ^ % + i - p i (5) 

# 

and, if equation (5) is substituted in equation (4), 

2 

% “ ^Si+i-*- ^ *1+1 " *1+1 " ^ p i + P i 

■ £ Ii+i + - ^) Ei+i + (* -£)**. (6) 

In addition, the following identity can he set up between the 
arbitrary radioactive powers of two elements: 

* 1-1 S ^~ *1 ( 7 ) 

JL 

with the following additional definitions: 



In matrix notation, equations (5) to (7) may be expressed as 
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The transmission matrix may "be abbreviated for conciseness of notation. 
Equation (8) then reads 



For shield element , similarly. 


Vl 

Xl*! 

P i-2 

The over-all effect of n elements is obtained by matrix multi - 
plication of transmission matrices: 
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where 


foi-l 

v p i, 




a - W 


rw 


P n+1 + P n 


where I + P is the power of radiation emitted from the shield 
n+i n 

to the surroundings and R^^ is the power of radiation emitted 

from the surroundings to the shield. From a practical aspect, 

Xj - + P n + R n+1 might "be the tolerance of a human being to radiation, 
and the background radiation in the vicinity of the Individual 

from sources other than 3L + 2 and ? n . 


If the source absorbs no radiation and reflects 


back 


to the shield, then 


II = *1 + Pi + Q 


where Q. is that portion of the power of radiation from the Bource 
intercepted by the area of the shield. If the source absorbs all 


back-scattered radiation, then 


is Q alone. The practical 


case lies somewhere between these extremes with the source absorbing 
some of the back-scattered power. 

Inversely, with and Ij known and and Ij +1 

desired, equations (5), (6), and (7) are rewritten as 


%+i - x ± - P i 


■ (>•!)», 


r i+l p 

TT 1 


I 
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where the definitions under equation (7) are now given with 
replacing P.j_-j_« In matrix notation. 




Por the entire shield. 


(16) 


\*1 


V 

V 

■ [a] 


0 


p l. 


(17) 


The over-all transmission matrix A is the same In this Inverse 

P 

circumstance as that previously given except that 1+1 replaces 

P. _ P 1 

I~1 .(provided that the order of B and I is reversed in the 

P i . 

column matrices), P-^ Is entered In the column matrix of the radia- 
tion entering the first shielding element, and P n Is removed from 
the column matrix of the radiation leaving the last element In the 
shield. The definitions given under equation (7) hold when P i+1 

replaces P^« 


Par a shield with no radioactivity, = 0 and a second-order 
transmission matrix may he set up from equation (8) : 



L ' . 

ti- ti 


*1+1 

bi H hj 


®i+l 


( 18 ) 
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The calculation of the over-all transmission matrix may "be performed 
"by matrix computation as indicated for the third -order matrix. In- 
asmuch as the reciprocal of this transmission matrix is the trans- 
mission matrix with the order of the rows and -the order of the 
columns reversed, the inverse problem may be set up with an identical 
transmission matrix provided only that the order of I ana r 
I s reversed: 


' B i + 1 ■ 


i 


V 

I i+1 


i 

•ril «H 
l-P 

1 

St?*. 


h 


If a complete shield is to be analyzed, the first step is to 
break down the thickness into a number of elements, the materials 
of which are known, for whose thicknesses transmission ana reflection 
data are available • The radiation from the source may be assumed 
to be monochromatic or to consist of a numb er of monochromatic 
components . This radiation will be assumed not to be degraded by 
passage through the shield but only to be absorbed and "back-scattered. 
I n asmuch as degradation of energy is more common for most applications 
of thick shielding than either upgrading of energy by nuclear reaction 
or tr an s m ission with energy unchanged, and inasmuch as neutron absorp- 
tion generally increases with decreasing energy (neglecting resonance) 
as does gamma-ray absorption (in the range of energy above the K- 
absorption edges and below the point at which pair formation causes 
the’ absorption to rise with increasing energy) , the assumption of 
no degradation of energy leads to the design of a safe shield. The 
shield may be less safe than indicated by analysis if the energy 
of radiation is above the energy at which- a mini mum exists in the 
absorption characteristic. If the power transmitted and "back- 
scattered by an element is known, application of the analysis will 
produce a very safe shield, in general, because the analysis not 
only neglects the effects of degradation of energy by scattering, 
but assumes a path for traverse of radiation normal to the shield, 
whereas actually the radiation traverses a longer path at some n»gi a 
to the , normal. 

Exangle 1. - Assume (l) that tj = 0.8 and r± = 0.1 
for all values of i at same energy for same material of a 
certain elementary thickness and that t and r are con- 
stants; (2) that the shield is 3 elements in thickness; (3) that 



12 


KACA-TH No. 1712 


the power of radioactivity from each element Is 0.0020 roentgen 
per hour and that the spectrum Is of the same type as that of the 
source; (4) that 1 background radiation is negligible; and (5) that 
human tolerance is 0.01 roentgen per hour. What size source can 
safely he shielded? 

These data are substituted in equation (ll) to yield 


-.2 


Pi-. " 


l 

0.1 

1 


1 


0.1 

1 

-^l 


0.8 

oTq 

0.8 


0.6 


0.8 0.8 

-n 


0.1 o 

B 0-01 

, 0.1 


OA n 

8- 

0.01 n 

0.1 

*i 



B 0.8 

0.8 


0.8 

0.8 

0.8 

0 


0 

0 

0 


0 


0 

1 













"9737 

20,169 

26 ,00l“ 


0.0100 


0.0089 




5120 

51,200 

5120 








20,169 

227,447 

82,863 


0.0000 

nr 

0.0072 




51,200 

512,000 

51,200 








0 

0 

0 


0.0020 


0 



0.0100 


0.0000 


0.0020 


Therefore 


1^ = 0.0089 roentgen per hour 
=- 0.0072 roentgen per hour 


and 


I x - El = 0.C017 roentgen per hour 

The source power Q equals if the source absorbs all 
radiation; Q = 1^ - ** ' the source reflects all radiation from 

the shield back to the shield. The result may consequently be 

expressed as / u 


0.0089 > Q ^ 0.0017 
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Example 2. - Assume (l) that t. = 0.8 and r.^ = 0.1 for all 

values of i at some energy for same material of a certain elementary 
thickness j (2) that the source power is 2 roentgens per hour; (3) that 
human tolerance is 0.0099 roentgen per hour; (4) that background radi- 
ation is negligible; and (5) that the material of the shield is not 
radioactive. How many elementary thicknesses are required to shield 
the source safely? 

Two alternatives may be set up, one in which the source absorbs 
all back-scattered radiation and one in which it absorbs none. 

When the source absorbs all back- scattered radiation. 




2 


~1 

0.1 " 


- 

0.8 

0.8 

■R 


0.1 

0 8- 0 * 01 



0.8 

°* 8 0.8 J 


0.0099" 


0.0099~ 


-H* 


0.0000 


0.0000_ 


By calculation. 


and 


M 26 

[*P 


0.0099 


1.6554 

p.oooo. 


_0.4859_ 




0.0099 


2.0085 

0.0000 


0.5896 


The result indicates that 27 elements of shielding are needed to 
shield a source with an intensity of 2 roentgens per hour if the 
source absorbs all back-scattered radiation from the shielding. 

When the source absorbs no back-scattered radiation. 
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2 + bJ 


1 

0^ 

l 

Him 
•1 • 

olo 

1 

z 

0.0099“ 


0.0099’ 

*L 

s 

0.1 

0.8-2^ 


0.0000 

= D*f 

0.0000 


0.8 

0,0 0^ 




By calculation, 



0.0099 


2.4369 


1.7215 + 0.7154 

„28 

M 

0.0000 

= 

0.7154 

CS 

0.7154 


[0.00991 


2.9567“ 


2.0887 + 0.86801 

El 29 

0.0000 

s 

0.8680 

= 

0.8680 


The result therefore Indicates that 29 elementary thicknesses of 
shielding are required to shield a source whose intensity is 
2 roentgens per hour if the source absorbs no radiation back- 
scattered from the shielding. 


Case II - Conversion to Beat 

The following hypothesis is made: 

Monochromatic radiation 
Radiation norma], to shield « 

Source of radiation outside shield 
Radiation absorbed converted to heat 

Method. - The conversion of radiation to heat is included; the 
amount of radiation absorbed is proportional to the amount incident 
upon the element. Radioactivity of the shielding is not considered. 
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The equations relating the powers of radiation incident on and 
scattered from the i**“ element are 


t i I i + Vi+l 


I i+1 

B i “ t i fi i+l + r i I i 


( 20 ) 


H i = Vi + Vi+i 


Solving for I*, and S» explicitly and setting in matrix form 

gives 

'l_ 

tj tj 


V 


H. 

s 

i 


h 



*1 


hi 

h 


_rj* 

hiri 

% 


I i+1 


1+1 

^+1 




L i+i 

^+1 

* 1+1 


( 21 ) 
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Ibr the complete shield. 


V 



^n+l 

% 

- 60 60 

• • • 

^l 

% 



_0 



I w-1 



ii 

r>i 

L_J 





0 



An alternate method of determining consists in determining 
only Ii and Rj by matrix methods and then commuting from the 
third of equations (20). 

Example 1. - Assume that an element of a uniform shield transmits 
0.8 of incident radiation, "back-scatters 0.1 of it, and absorbs the 
rest; the tolerance level is 0.0099 roentgen per hour. If the source 
power is 0.1626 roentgen per hour, how much heat does a ld-element 
shield absorb? How much heat does the source absorb? 

Die values of the power coefficients are 


t « 0.8 
r = 0.1 

h = 1 - t - r = 1 - 0.8 - 0.1 = 0.1 
Application of equation (22) yields the following data: 


V 


0.0099 

% 

- m u 

0.0000 

% 


_ 0 


10.1626 

1 0.0476 (22a) 

[0.0202 

Successive applications of equation (21) yield the following table: 
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Element 

Heat absorbed, 
(equivalent 
roentgens/hr) 

14 

0.0012 

13 

.0017 

12 

.0021 

11 

.0027 - 

10 

.0034 

9 

.0041 

8 

.0051 

.7 

.0062 

6 

.0078 

5 

.0093 

4 

.0113 

3 

.0137 

2 

.0166 

1 

.0202 

0.1052 total heat absorbed in 

shield 


The amount of radiation absorbed as heat In the source 1 b the 
amount hack-scattered from the first element of the shield B^; 

= 0.0476 roentgen per hour hy calculation from equation (22a). 


Case HI - Polychromatic Radiation and Degradation of 

Energy 

The following hypothesis is made: 

Polychromatic radiation 

Radiation normal to shield 

Source of radiation outside shield 

Radiation absorbed converted to heat and mass 

Radiation degraded or upgraded In energy 

Method. - The radiation incident on and hack-scattered from an 
element of shielding Is assumed to consist of a spectrum of radiation 
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The power coefficients t^ and take on different values for 

each portion of the spectrum incident on the element of the shield. 

For example, at an energy of 2.0 Mev, t^ mi g h t he 0.82 for 

whereas, at 1.8 Ifev, tj might he 0.81. In the same way, the pro- 
portions of radiation converted to heat and mass depend on the energy 
of radiation incident on the element. At 0.9 Mev, for example, 

= 0 for gamma rays because no pair formation occurs. 

The power coefficients tj, rj , hj, and must consequently, 

he regarded as operators that convert one spectrum of radiation to 
another, such that the coefficients at each portion of the spectrum 
alter the corresponding intensities differently. 

Some of the energy also must he degraded (neglecting nuclear 
reactions) after encountering an element of the shield. Thus, if 
a monochromatic "beam of gamma radiation of 1 Mev encounters an 
element of a shield, much of the l-Mav radiation passes through 
undisturbed in energy, hut some radiation is degraded to the range 
of, perhaps, 0.9 to 1 Ifev, some to the range 0.8 to 0.9 Jfev, and 
so forth. For these calculations, it is convenient to regard an 
incident spectrum of radiation as consisting of a number of dis- 
crete energy ranges, each of a certain intensity. The power 
operators then convert each range to a number of other ranges of 
energy besides permitting some of the origi n al range to pass 
through. 

These operators may he set up as matrices. Let the incident 
spectrum of radiation he represented as 


F(x) - a^l + a 2 x 2 + 


+ + 


• • • 


+ 
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where x^ represents unit power at energy i. Then the array of 
intensity coefficients nay "be considered to he a matrix: 



The transmission matrix may he written correspondingly in expanded 
matrix form: 


t = 


*11 *12 * 

• 

• 

*11 * 

• • 


*21 *22 * 

• 

• 

*21 * 

• • 

^ 2n 


*11 t i2 




*nl *n2 


-*ni 


nn 


(24) 


where tj i is that part of the operation that converts incident 
energy Xj to emergent energy x.. The operation tF may he 
performed hy matrix multiplication. 


An examination of the matrix t reveals that the terms above 
the principal diagonal represent upgrading of energy as hy nuclear 
reaction; the terms helow the main diagonal represent degrading of 
energy; and the terms on the main diagonal represent transmission 
of radiation of unchanged energy. _A common case is that for whi ch 
no upgrading occurs. The matrix t tak es on trian gular form for 
that case. ■ The operator h contains terms only on the pr inc i pal 
diagonal. 
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Example. - Let a hypothetical spectrum of radiation of the 
following energy distribution he incident on a shield element: 


Range of 
incident 
energy 
(Mav) 

Relative power 
of incident 
radiation 

1.0-0. 9 

10 

.9- .8 

9 

.8- .7 

8 

.7- .6 

8 

.6- .5 

7 

.5- .4 

6 

.4- .3 

6 

.3- .2 

7 

.2- .1 

8 

.1- .01 

9 

Total 

78 


This distribution appears in matrix form from equation (23) : 


F 


10 

9- 

8 

8 

7 

6 

6 

7 

8 
9 


Also given are the transmission, reflection, and absorption character- 
istics of a hypothetical shield element as shown in the following table: 
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22 


MCA TN HO. 1712 


The problem Is to find the intensity distribution of idle radiation 
transmitted and back-scattered, and to determine the power of radiation 
absorbed as heat. 

From the data given, 

— 1 T~ 

10 

9 

8 

8 

7 

6 

6 

7 

8 
•9 


tF = 


0.70 0 0 

0 

0 0 

0 

0 

0 


0 

0.01 0.65 0 

0 

0 0 

0 

0 

0 

• 

0 

0.01 0.01 0.60 

0 

0 0 

0 

0 

0 


0 

0.01 0.01 0.01 

0.50 

0 0 

0 

0 

0 


0 

0.01 0.01 0.01 

0.01 

0.40 0 

0 

0 

0 


0 

0.02 0.02 0.02 

0.02 

0.03 0.30 

0 

0 

0 


0 

0.01 0.02 0.02 

0.02 

0.02 0.03 0.20 

0 

0 


0 

0.01 0.01 0.02 

0.02 

0.02 0.03 

0.04 

0.15 

0 


0 

0.01 0.01 0.01 

0.02 

0.02 0.02 

0.03 

0.05 

0. 

.15 

0 

_0.01 0.01 0.01 

0.01 

0.01 0.02 

0.03 

0.05 

0.05 

0.10 

*7.0040 4-0 4- 

0 + 

0 4- 0 + 

0 + 

0 + 

0 

+ 

0 

0.1045.85+0 4- 

0 + 

0 + 0 4- 

0 + 

0 + 

0 

+ 

0 

0.1040.09+4.80+0 + 

0 + 0 + 

o + 

0 + 

0 

+ 

0 

0.1040.09+0.0844.00+0 + 0 + 

0 + 

0 + 

0 

+ 

0 

0el0-*0#09+0#08+0. 08+2 *80+0 + 

0 +■ 

0 +. 

0 

+ 

0 

0.2040.18+0.1640.16 40. 21+1. 80+O - + 

0 + 

0 

+ 

0 

0.1040.1840.16+0.1640.1440.1841.2040 +- 

0 

+ 

0 

0.10+0.09 40.1640.1640.1440.1840.24+1.05+0 

+ 

0 

0.1040.0940.0840.1640.1440.1240.1840.35+1.20+0 

0.1040.0940.0840.0840.0740.12+0.1840.35+0.40+0.90 


The operations rF and hF may be performed in a similar manner. 
The results are summarized in the following table: 
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mlttafl.. back- acattarad . and abaortiea 



.00 1.60 


0.90 

1.80 


8.71 

1.41 


8.18 

1.97 


8.18 

3.15 


8.48 

3.38 
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Equations may "be written for the disposition of radiation inci- 
dent on a shield element, as in previous cases. Die figure illus- 
trating case HI gives 


h.+i = + r i B i+il 


% = %%+! + r i I i 
% = ^L 1 ! + “A+l 
% = + SiEi+l^ 


(25) 


__ _In the interpretation of these equations, the fact that t, r, 
m, h, 3j_, Xj_ +1 , Ej, Mp and are matrices must he 

kept in mind. 3n general, two matrices are noncommutative in multi - 
plication; t^ ^ %^i-i3i except for certain types of matrix 

if tj- and t^_ 2 . are different. Practically, this nancamrautativity 
means that an element of carton in front of an element of ledd, in 
general, would act differently on radiation from an element of lead 
in front of an element of carbon. 


Before transposition of terms in equations (25), an examination 
into the existence of the reciprocal matrix «na into the uniqueness 
of solution for the inverse operation when t is a singular matrix 
must he made. Die reciprocal matrix exists if, and only if, the 
determinant of the original matrix is nonzero. Die conditions of 
the problem state that no energy level is completely extinguished 
and consequently the diagonal elements of "the original matrix are 
nonzero. A frequently occurring transmission matrix is -that cor- 
responding to the case of no upgrading of energy; the elements of 
the matrix above the principal diagonal are all zero so that the 
matrix is triangular. Inasmuch as a triangular matrix with all 
nonzero elements on the diagonal is never singular, the reciprocal 
matrix in this instance always exists. 

A transmission matrix may he singular, then, only when upgrading 
of energy takes place in the transmission process. Diere are two 
practical methods by which a solution may he obtained when the trans- 
mission matrix is singular. Die first method (reference 3) yields 
a solution in terms of a number of arbitrary parameters equal to the 
degeneracy of. the matrix. Diese parameters, together with the boundary 
condition imposed by a knowledge of ihe spectrum incident on the 
first element, may yield a unique solution to Ihe problem. Die 
second method consists of Ihe arbitrar y adjustment of the elements 
of the transmission matrix so that it is no longer singular. Die 
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number of elements to be adjusted is, in general, equal to the degen- 
eracy of the transmission matrix. QMs distortion of data is insig- 
nificant if the adjustments are smaller than the experimental errors 
in determining the data. 

Multiplication from the left of equations (25) by t. (the 
reciprocal of the matrix tj) and transposing yields 


= % 1 -^i+l " % r i B i+l 

% - Vf 1 ^i+i + - ? i V 1 V Vi 

% = SlV 1 *1+1 + (^ - % tj" 1 ? ± ) 1 (26) 

% = ^i^i I i+1 + 0»i - *i t ± r ± ) %+! 

If equations (26) are rewritten in matrix notation and if the 

bars over the symbols for the operators are dropped for convenience 
of notation. 


% . r i 0 0 *1+1 

Vi" 1 % - Vi" 1 0 0 %+i 


“-V 1 r i 


Vi " 1 “i - Vi " 1 r i 0 0 %+i 

*1% 1 . hj: - hjtj 1 rj 0 0 H^ +1 
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and 


X 

H W • • • 


I n*l 

^n+1 

-H 

W 

^n+l 


% 



0 


0 

% 



0 


0 


As usual, the sum of the elements In the matrices 


(28) 


i=n 

% 

1=1 

is the total heat absorbed "by the shield* 

It may he noted that a spectrum can he chosen for that 

could not possibly have been emitted by a given shield upon which 
a spectrum 1^ was incident. The following exangale will be examined 
for illustration: Let = 0. let X^ + i and Ij_ each consist 

of two energy ranges* Then 



let 


t i I i " I i+1 


t 


0.4 0 
0.4 0.9 


By substitution in equation (28a) 


0.4 


0.4 


0.9 


Z n 


12 1 


°- 4x H 


O.Az^ + 0.9 z 12 


^i+1 


(28a) 
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Obviously, I^ + 2 _ can never be chosen to be a spectrum whose upper- 
energy region is of higher power than its lower-energy region, because 
such a choice leads to the absurdity of a negative power for some 
range of energy in 1^. The problem of shielding polychromatic radia- 
tion must be solved by starting from only the incident radiation 
instead of from the emitted radiation in order to avoid such absurd 
results. 

let equation (28) be rewritten: 


h 


% 

U 12 ^ 


W 


*1 

- 

"a 

U 22 U 23 



(29) 

<a 
— — 


v 

U 32- U 33_ 


_0 _ 



where the over-all transmission matrix A has been partitioned into 
submatrices U. .. Then 


or 


*1 = ¥ml + *12*n + l . 
E 1 " E 21 I n+-l + U 22 B n*l 


I n+1 = U 11 1 X 1 “ ^ll" 1 U 12 B ih-1 


®1 “ ^21^11 1 \ + ^22 " U 21 U 11 J ‘ U 12^&H 


-1 


(30) 


which is the required solution provided that the inverse matrices exist, 
inasmuch as all terms on the right side of equation (30) are positive 
or zero under the conditions of the problem. When the inverse matrices 
do not exist because the original matrices were singular, the methods 
of solution in terms of parameters or adjustment of the original 
matrices may be used. 


Case 17 - Gamma Plus Neutron Radiation 
The following hypothesis is made: 

Polychromatic gamma plus neutron radiation 



28 


HACA 037 No. 1712 


Radiation normal to shield 

Source of radiation outside shield 

Radiation absorbed converted to energy other than gamma 
rays and neutrons 

Gamma radiation gives rise to neutrons and neutrons give 
rise to gamma rays 

Radiation degraded or upgraded in energy 

Mathod . - Incident gamma radiation 1^ and neutron radiation 

3? are considered to he partly transmitted (with and without energy 
change) , partly hack-scattered, partly converted to mass and heat, 
and partly converted one to the other (7 to n and n to 7 ). In 
order to minimi ze the amount of notation, both (radiation power 
converted to mass other than neutrons) and (radiation power 

converted to heat) will he combined into one term C. (converted 
power) . 



% + = Cj 


3he following equations may then he set up (the bars over the 

operators t, r, and c will he dropped for convenience of 
notation in this case) : 


9 83 
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4* - 44 * <T=4 ♦ 

(31) 

» 

n nn yn y nn yn y 

*1+1 ° + % H + r i B i+l + r i B l +1 

(32) 

P 7 = +V 4 . + n V 4 . -V 4 - ^ n 7 T n 

E i t i R i+l + *i ^+1 + 'i 3 ! + r i ^ 

(33) 

n nn yn 7 nn yn y 

% " t i R i+l + t i %+! + 'i 3 -! + r i h. 

(34) 

°i - «1 <4 ♦ ®i+l) +<=?(!? + E?+l) 

(35) 


Equations (28) and (29) may he solved simultaneously for 1^ ana 
I^. These values may he substituted in equations (30) to (32) to 
yield values of E^, E^, and in terms of only l£ +1 , 

■^i+l^ ^i+i* and E^+ 2 * She resulting matrix has the farm: 


A 


Z 11 z 12 z 13 z 14 z 15 

*1 «' 


A 


Z 21 Z 22 Z 23 Z 24 *25 

ijl 

i+1 


E i 

= 

*31 Z 32 z 33 *34 *35 

R i+1 


* 


z z z z z 
41 42 43 44 45 

Cl 


_ C i_ 


Z 51 *52 z 53 z 54 z 55 

— Ml 

°i+l 


A 


Cl 


4n 

3 i 


*1 


4n 

*1 

= 

hi H • • • &>] 

Ci 

-M 

®Li 



Ci 


Cl 

_pi_ 


_ 0 


. 0 _ 


which is the required solution. 


(36) 


( 37 ) 
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The calculations here are lengthy "because of frequent occurrence 
of inverse operators. Machine methods are necessary for handling a 
case of this complexity. 


Case V - Eadiation Not Normal to Shield 
The following hypothesis is made: 

Polychromatic radiation 

Radiation not normal to shield (spherical shield) 

Source of radiation outside shield 
Radiation absorption neglected 
Radiation degraded in energy 
Inverse -square decrease in power neglected 

A spherical shield was chosen because reasons of symmetry all ow 
neglect of the displacement of radiation from the normal to the shield 
that is caused by scattering at an angle other than 180°. In order 
to handle directional scattering, the incident radiation may be 
assumed to arrive at a unit area of the shield Sj in a number of 
conical zones identified as a and P when two zones are chosen. 

An equal number of zones of emergence are chosen and designated by 
primes. 

a 


P 

a 


The equation for scattered radiation may consequently be set 

up: 



I i+l,a' = \a ,I i,a + ^pa ,:E i,P + r p'a' R i+l,p* + r a'a ,B i+l,a* 

Expressions for Ij + x,| 3 «j 2114 "be similarly set up 

in matrix form. The coefficients are not ail" independent (neglecting 
Bragg reflections) : 
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.■kpa* = ^3’ 
r Pa = r p«a* 

By consideration of the difference in path length in traversing the 
element at angles greater than 90°, teg, and can he related. 

The coefficients for the path lengths that are not; integral multiples 
of the thi ckn ess of the standard element can he obtained by deter- 
mining a fractional power of the transmission matrix. The question 
arises as to whether this fractional, power of a matrix is unique 
for shield in g applications where all elements of matrices must he 
positive and real, and for special matrices such as tr iangular matrices 
that result when energy is not upgraded. These and other questions 
" that must he considered when analytically determining the transmission 
matrix for an element that is not an integral multiple of the st andar d 
thickness will not he considered because it was assumed that the data 
could he experimentally obtained. 

A question that must he considered is whether degeneracy exisus 
to any greater extent in transmission matrices of the t a p, and r a g 
types than in the t^, type, and what are the consequences of this 
greater degeneracy, if any. It is physically apparent that inasmuch 
as radiation in the incident and emergent beams are no longer in the 
same direction far the t^p, type of transmission, no condition on 
the problem prevents complete extinguishing. 3h two cases, however, 
this extinguishing will not occur: (1) When adequate amount of 

unmodified scattering is obtained; and (2) when modified scattering 
does not degrade energy out of initial energy range (with upgrading 
of energy neglected) • For these cases, the transmission matrix is 
never singular. In general, however, more sin gular transmission 

matrices of the t^, type may he expected than of the t , type. 


DISCUSSION 

A l imi ted number of cases have been presented to indicate the 
wide variety of problems that may be handled by matrix methods. 

Die inverse-square decrease in intensity with distance from the 
source has been neglected. This decrease can be taken into account 
by multiplying the radiation power at any station in the shield by 
the ratio of the area of the first element to the area of the ele- • 
ment at the given station. Ohus, if I^ and B-j_ were determined 
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at the surface of a sphere 2 feet in diameter and 3^ and at 
a sphere 6 feet in diameter, the calculated values of I n and 
would "be multiplied by l/9 to include the inverse-square decrease. 

Data on the directional distribution of intensity and energy at 
high energy are required for gamma rays. Ifeny available theoretical 
expressions for scattering and absorption may be extrapolated to 
* higher energies to provide -these data. Equations such as the Klein- 
Uishina equations (reference 4) appear to be adequate, at least in 
the range of energy that has been experimentally checked. Experi- 
mental data exist in the form of tables, particularly for wavelengths 
above 0.005 A. One of the most comprehensive compilations of mass- 
absorption coefficients for high-energy radiations ■ may be found in 
reference 5. Such mass-absorption coefficients yield information on# 
the amounts of radiation transmitted in a primary beam through a thin 
absorber for many elements and a number of wavelengths above 0.005 A. 
The radiation considered absorbed Includes not only such radiation 
as may be converted to heat, but also radiation diffracted by the. 
crystal lattice of the absorber, radiation deviated from the primary 
beam by Compton scattering, and same radiation participating in pair 
formation (a factor almost negligible below the maximum energy of 
2.5 Mev tabulated). Combinations of experimental and theoretical 
data yield truer results than either alone. 

Data are required only for scattering and absorption by chemical 
elements and not by compounds and mixtures because Interaction between 
nuclear radiation and solids is negligibly affected by binding forces 
between atoms or by the physical state of the element. The elementary 
thickness of absorber may consequently be regarded as that thickness 
that contains a weight of element per unit area that is regarded as 
standard. For example, if 2 grams per square centimeter of element 
A and 1 gram per square centimeter of element B are standard where 
the atomic weight of B is, say, twice that of A, then the com- 
. pound AB in thickness g grams per square centimeter may be 
regarded as containing alternately standard thickness of A and B 
for a total thickness of g/3 elements of thickness each of A and B. 

The inherent advantage of requiring data for only one thick- 
ness of absorber is tempered by the fact that some thicknesses may 
be inadequate for accurately determining variation of spectral 
distribution of energy with angle of scattering. Different thick- 
nesses of absorber may be required to obtain experimental data in 
same instances. 

It should be noted that elements of low atomic number do not 
necessarily possess poorer mass-absorption coefficients than elements 
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of high atomic number for X-rays and gamma rays despite the fact 
that the low-number elements do absorb most poorly at wavelengths 
of 0.02 A and below (neglecting absorption edges). Hydrogen is a 
notable exception of this sort Inasmuch as the data of reference 5 
show that hydrogen possesses a mass-absorption coefficient 64 percent 
greater than that of lead at a radiation of 0.010- A wavelength, and 
83 percent greater at 0.005 A. At 0.005 A, the coefficient of carbon 
is only 10 percent lower than that of lead. At 0.005 A, a long- 
chain hydrocarbon that might be regarded as possessing essentially 
twice as many hydrogen atoms as carbon atoms would possess two parts 
by weight that absorb 83 percent more than lead, and 12 parts by 
weight that absorb 10 percent less than lead — a combination that 
absorbs gamma rays as well as lead, weight for weight? 

Machine methods for solution of matrix multiplications are 
certainly necessary for comprehensive solutions. If, for example, 
both neutron and gamma radiation are considered, if heat and mass 
conversion are desired, and if five zones of scattering are 
chosen, a 3 Oth- order square matrix is set up, each term of which 
is a matrix of order equal to the number of energy zones into 
which the spectrum is divided. For a 5-foot thickness of shield- 
ing, 60 matrix multiplications are required if the element is 
1 inch in thickness; 600 if the element is 0.1 inch in thickness. 

The type of composite shield that should be used would be 
determined by trlal-and- error calculations. A single material 
might be used for a first computation. TJjpon analysis of that 
hypothesis, observation might reveal that the spectral distribution 
of neutrons and gamma radiation is such at any station in the 
shield that the succeeding elements might well be another material 
better fitted for absorbing or degrading energy. The calculation 
would then be repeated and the result reanalyzed. Such calcula- 
tions, although lengthy for accurate results, might be more eco- 
nomical than experimental verification of the adequacy of complete 
shields. 

A number of simplifications to the exact treatment of the 
shielding problem can be made by making approximations on some of 
the shielding processes. The following table lists same changes 
that readily suggest themselves in order to reduce the amount of 
confutation (all the suggested approximations yield an oversafe 
shield, that is, thicker or heavier than required): 
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Exact mechanisms 

Approximating assumptions 

Radiation is scattered con- 
tinuously at angles to pri- 
mary direction of incident 
radiation 

That small number of finite 
zones exist; that radiation 
may be allocated only to zones 
having equal or lesser path 
length through shield element 
than in true case; that all 
radiation scattered passes 
normally through element. 

Spectrum of radiation is 
scattered 

That all radiation is mono- 
chromatic, if possible of 
that energy at •which minimum 
absorption occurs. 

Radiation is degraded or 
acted upon differently 
according to energy 

That radiation is degraded or 
acted upon by, constant factor 
independent of energy. 


CONCLUSION 

An analysis of the shielding of radiation indicates that 
matrix methods appear to he immediately applicable for solution 
of shielding problems, provided that adequate experimental data 
on scattering and absorption for standard thickness of absorbers 
are available. 


Plight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, June 7, 1948. 
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